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Abstract 

The dynamic critical behavior of the two-rephca cluster algorithm is studied. Several versions 
of the algorithm are applied to the two-dimensional, square lattice Ising model with a staggered 
field. The dynamic exponent for the full algorithm is found to be less than 0.4. It is found that 
odd translations of one replica with respect to the other together with global flips are essential for 
obtaining a small value of the dynamic exponent. 
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I. INTRODUCTION 



The Swendsen-Wang (SW) algorithm and related cluster methodd^'lffli'l have greatly 
improved the efficiency of simulating the critical region of a variety of spin models. The 
original SW algorithm can be modified to work for spin systems with internal symmetry 
breaking fieldJ. Spin models of this kind include the Ising antiferromagnet in a uniform 
field, the random field Ising model and lattice gas models of adsorption in porous medial. 
The modification proposed in Ref. ^ is to assign Boltzmann weights depending on the net 
field acting on the cluster to decide whether the cluster should be flipped. Unfortunately, the 
modified SW algorithm is not efficient. The problem is that large clusters of spins usually 
have a large net field acting on them and are prevented from flipping by these fields. An 
algorithm for Ising systems with fields that avoids this problem was introduced by Redner, 
Machta, and Chayeffl. In this two-replica cluster algorithm large clusters are constructed 
from two replicas of the same system and have no net field acting on them so that they 
may be freely flipped. The two-replica cluster algorithm has been applied to study the 
phase transition of benzene adsorbed in zeolites! and is more efficient than the conventional 
Metropolis algorithm for locating and simulating the critical point and the phase coexistence 
line. Combined with the replica exchange method of Swendsen and WangEl, the two-replica 
method has been applied to the random field Ising modelS The two-replica method is 
closely related to the geometric cluster Monte Carlo methodBSEl. 

In this paper, we report on a detailed investigation of the dynamics of the two-replica 
cluster (TRC) algorithm as applied to the two-dimensional Ising ferromagnetic in a stag- 
gered field (equivalently, the Ising antiferromagnet in a uniform field). The TRC algorithm 
introduced in Ref. |10| has two components that are not required for detailed balance and er- 
godicity. We studied the contribution to the performance of the algorithm of these optional 
components. We find that the complete TRC algorithm has a very small dynamic expo- 
nent z < 0.4. However, we also find that this small value of z requires one of the optional 
components and that this component depends on a special symmetry of Ising model in a 
staggered field. This observation leads to the question of whether cluster methods exist for 
efficiently simulating more general Ising models with fields. We investigated other optional 
components for the algorithm but these do not lead to acceleration when fields are present. 

This paper is organized as follows. In Sec. II we introduce the Ising model in a staggered 
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field and describe the TRC algorithm. In Sec. |T| we define the quantities to be measured 
and how errors are computed. In Sec. fV| we present the results. The paper closes in Sec. [V] 
with a discussion. 



II. THE MODEL AND TWO-REPLICA ALGORITHM 

A. Ising Model in a Staggered Field 

The Hamiltonian for the Ising model in a staggered field is 

(3H[a] = -K u,G, - H,a, (1) 

<i,j> i 

where the spin variables, o"j take the values ±1. i^' is the coupling strength and Hi is the 
magnetic field at site i. The summation in the first term of Eq. (|1]) is over nearest neighbors 
on an L X L square lattice with periodic boundary conditions and L even. The second 
summation is over the sites of the lattice. The staggered field is obtained by setting Hi = H 
if i is in the even sublattice and Hi = —H if i is in the odd sublattice. The staggered 
field breaks the up-down symmetry((Tj ^ — cTj) of the zero field Ising model, however two 
symmetries remain. The Hamiltonian is invariant under even translations: 

ai+ro ^ ^ (2) 

with To any vector in the even sublattice. The Hamiltonian is also invariant under odd 
translations together with a global fiip: 

cTi+n < o"j for all i (3) 

with ri any vector in the odd sublattice. 

Figure ^shows the line of critical points, Kc{H) for this model. We carried out simulations 



at three points on the critical line taken from the high precision results of Ref. ^ 

KciO) = 0.4406867952 

Kci2) = 0.7039642053 
K^{4:) = 1.1717153065 
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The basic idea of the two-rephca cluster algorithm is to simultaneously simulate two 
independent Ising systems, a and r, on the same lattice and in the same field. Clusters of 
pairs of spins in this two-replica system are identified and fiipped. In order to construct 
clusters, auxilliary bond variables are introduced. The bond variables {rjij} are defined for 
each bond < i,j > and take values and 1. We say that < i, j > is occupied if rjij = 1. A 
bond < i,j > is satisfied if (Xj = aj and Xj = tj. Only satisfied bonds may be occupied. 

The two-replica algorithm simulates a joint distribution of the Edwards-Sokal0 type for 
{cTj} and {tj}, and {rjij}. The statistical weight X[a,T,ri] for the joint distribution is 

X[a,r,r7] = e-^['^'^A[a,r,r/]5pM (4) 

where 

<i,j> i 

B is the standard Bernoulli factor, 

Bp[iT\=p\'^\l-pY'^-\^\ (6) 

\rj\ = ^ {< i, j > \riij = 1} is the number of occupied bonds and Nb is the total number of 
bonds of the lattice. The A factor enforces the rule that only satisfied bonds are occupied: 
if for every bond < i,j > such that rjij = 1 the spins agree in both replicas (cTj^ = aij 
and Tj = Tj) then A[a,T,ri] = 1; otherwise A[a,T,ri] = 0. It is straightforward to show 
that integrating X[cr, r, 77] over the bond variables, rj yields the statistical weight for two 
independent Ising model in the same field, 

^-m[^]-m[r] ^ ^^^gt X[a, r, ri] (7) 
if the identification is made that p = 1 — e~^^ . 

B. Two-Replica Cluster Algorithms 

The idea of the two-replica cluster algorithm is to carry out moves on the spin and bond 
variables that satisfy detailed balance and are ergodic with respect to the joint distribution 
of Eq. (^. The occupied bonds 77 define connected clusters of sites. We call site i an active 
site if (jj 7^ Tj and clusters are composed either entirely of active or inactive sites. If a cluster 
of active sites is fiipped so that a < a and r < r the factor G is unchanged. 
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A single Monte Carlo sweep of the TRC algorithm is composed of the following three 
steps: 

1. Occupy satisfied bond connecting active sites with probability p = 1 — e~^^ . Identify 
clusters of active sites connected by occupied bond (including single active sites). For 
each cluster /c, randomly and independently assign a spin value Sk = ±1- If site i is in 

cluster k then the new spin values are cTj Sk and Tj < s^- In this way all active 

sites are updated. 

2. Update each replica separately with one sweep of the Metropolis algorithm. 

3. Translate the r replica by a random amount relative to the a replica. If the translation 
is by an odd vector, all r spins are flipped. 

Step 1 of the TRC is similar to a sweep of the SW algorithm except that clusters are 
grown in a two-replica system rather than in a single replica and only active clusters are 
flipped. Note also that the bond occupation probability is p = 1 — e~^^ for the TRC 
algorithm and p = 1 — e~^^ for the SW algorithm. It is straightforward to show that Step 
1 of the TRC algorithm satisfies detailed balance with respect to the joint distribution Eq. 
(|). Since only active sites participate in Step 1 of the algorithm, the Metropolis sweep, 
Step 2, is required for ergodicity. Step 3 contains the optional components of the algorithm: 
an even translation or an odd translation plus flip of one replica relative to the other. These 
moves are justified by the symmetries of the Ising model in a staggered field stated in Eqs. 
(0) and d^). When we refer to the TRC algorithm without further specification, we mean 
the algorithm described by the Steps 1-3 above. In the foregoing we also study the TRC 
with only even translations or with only odd translations. 

In the TRC algorithm we flip only active clusters but it is also possible to flip inactive 
clusters if a weight factor associated with the change in G is used. We call a flip of an active 
cluster to an active cluster (h — to — h or — h to H — ) an active flip. The TRC algorithm 
with inactive flips is obtained by replacing Step 1 with the following: 

1'. Occupy satisfied bonds with probability p = 1 — e~'^^ . Identify clusters connected 
by occupied bonds (including single sites). For each cluster k, taken one at a time, 
randomly propose two new spin values values, = ±1 and t^ = ±1 for the cr and r 



5 



spins respectively. Compute 5G, the change in G that would occur if the spins in the 
k*^ cluster are changed to the proposed values leaving spins in other clusters fixed. If 
5G < accept the proposed spin values (set (jj Sk and Xj ^ tk for all sites i in 
cluster fc), otherwise, if 5G > accept the proposed spin values with probability e"*^*^. 

Step 1' is by itself ergodic however it may be useful to add Metropolis sweeps and trans- 
lations. 



III. METHODS 

We measured three observables using the TRC algorithm: the absolute value of the 
magnetization of a single replica, m; the energy of a single replica, S; and the absolute value 
of the net staggered magnetization for both replicas, s, where the definition of s is 

^=I(E - E + (8) 

i£odd ideven 

Note that the staggered magnetization is conserved by all components of the TRC algo- 
rithm except Metropolis sweeps and inactive fiips. For each of these observables we com- 
puted expectation values of the integrated autocorrelation time. Tint and the exponential 
autocorrelation time, r^xp- From Tint, we estimated the dynamic exponent z. 
The autocorrelation function for 0, T^^{t) is given by, 

P = „,„ £K.Wf)-^)W'' + ')-^). (9) 

The integrated autocorrelation time for observable is defined by 

^ = ^+1™ Er<^0W (10) 

and the exponential autocorrelation time for an observable is defined bylll 

Texp,4> = lim - — tJ^ttt- (11) 

t-oo log r^<^(t) 

In practice the limits in Eqs. (P), (0) and (0) must be evaluated at finite values. The 
length of the Monte Carlo runs determine / and are discussed below. Following Ref. 19, we 
define 

r^nt,<p=l + j2^Hit) (12) 
^ t=l 



6 



and choose the cutoff t* to be the smallest integer such that t* > KTint,<j,, where k = 6. We 
used the least-squares method to fit log r0(^(t/rj„i_<^) as a function of t to obtain the ratio of 
rint,^/Texp,<t, and chose a cut-off at t/rint,^ = 5. 

We used the blocking methodB'El to estimate errors. The whole sample of n MC mea- 
surements was divided into m blocks of equal length / = n/m. For each block i and each 
measured quantity A, we computed the mean Ai . Our estimates of A and its error 6 A are 
obtained from: 

A=-Y.A, (13) 

1 m 

SA^= , ^, yiA-A,f (14) 

In our simulations, we divided the whole sample into m blocks where m is between 10 and 
30. 

For the data collected using the TRC algorithm, each block has a length I > 10^rj„(. 
For the data collected using modifications of the TRC algorithm, each block has a length 
/ > lO'^Tint- Data were collected for H = 0, 2 and 4 and for size L in the range 16 to 256. 

IV. RESULTS 

A. Integrated Autocorrelation Time 

Table | gives the integrated autocorrelation time using the TRC algorithm for the mag- 
netization, energy and staggered magnetization. Table | is comparable to the Table in Ref. 



T0| but the present numbers are systematically larger, especially at the larger system sizes. 
This discrepancy may be due to the sliding cut-off t* used here instead of a fixed cut-off at 
200 employed in Ref. |Tg. 

Table |l| gives the integrated autocorrelation times for magnetization using the TRC 
with only even or only odd translations. The comparison of TRC algorithm with only even 
translations and with only odd translations in Tables || shows that odd translations together 
with global flips of one replica relative to another are essential to achieve small and slowly 
growing autocorrelation times when the staggered field is present. 



Table 111 shows the magnetization autocorrelation times using different algorithms for 



system size L = 80. The Swendsen-Wang (SW) algorithm has the smallest Tint,m in the 
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absence of fields. However, when fields are present and the SW algorithm is then modified 
according to the method of Ref. ^ the performance is worse even than that of the Metropolis 
algorithm. The slow equilibration of the SW algorithm in the presence of the staggered field 
is due to small acceptance probabilities for flipping large clusters. On the other hand, the 
presence of staggered fields does not significantly change the performance the two-replica 
algorithm so long as odd translations are present. Inactive flips are helpful when there is 
no staggered field but when the staggered field is turned on, the autocorrelation time is not 
substantially improved by inactive flips. The explanation for the ineffectiveness of inactive 
flips when the staggered field is present is that the probability of accepting an inactive flip 
is small. For example, this probability is 1.4% for L = 80 and H = 4. 

The CPU time per spin on a Pentium III 450 MHz machine was also measured for the 



various algorithms and is listed in Table |T| for L = 80 . By considering a range of system 
sizes we found that the CPU time for one MC sweep of the TRC algorithm increases nearly 
linearly with the number of spins. The TRC algorithm is a factor of 3 slower than the 
Metropolis algorithm but this difference is more than compensated for by system size 80 
by the much faster equilibration of the TRC algorithm. Even without odd translations, the 
TRC algorithm outperforms Metropolis for size 80. 



B. Exponential Autocorrelation Time 

The ratio of the integrated to exponential autocorrelation times was found to be nearly 
independent of system size over the range L = 16 to L = 256. We found that over this size 
range Tint,m/rexp,m varied from 0.448 ± 0.008 to 0.425 ± 0.008 for H = 0; from 0.44 ± 0.01 
to 0.43 ± 0.01 for H = 2; and from 0.448 ± 0.009 to 0.409 ± 0.009 for H = A. The ratio 
Tint,s/'Texp,s is also nearly independent of L and H and is about 0.45. The ratio Tint,e / 'Texp,e is 
nearly independent of L but decreases slowly with H ranging from 0.29 to 0.25 as H ranges 
from to 4. The almost constant Tint,4>lTexp,4> for different sizes suggests that the integrated 
and exponential autocorrelation times are governed by the same dynamic critical exponent. 
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C. Dynamic Exponent 



Figures |] and ^ show the magnetization integrated autocorrelation time for the TRC 
plotted on log-log and log-linear scales, respectively. Figures ^ and | show the energy inte- 
grated autocorrelation time for the TRC plotted on log-log and log-linear scales, respectively. 
Figures ^ and |7| show the staggered magnetization integrated autocorrelation time for the 
TRC plotted on log-log and log-linear scales, respectively. 

For the whole range of L, logarithmic growth appears to give a somewhat better fit than 
a simple power law, particularly for the magnetization. Therefore, our results are consistent 
with z = for the TRC algorithm. Under the assumption that the dynamic exponent is not 
zero, we also carried out weighted least-squares fits to the form AL^ and varied L^m, the 
minimum system size included in the fit. Figures ||, ^ and |10| show the dynamic exponent 
z for the magnetization, energy and staggered magnetization, respectively, as a function 
of Lmin using the TRC algorithm. Figures and 0, show the dynamic exponent as a 
function of Lmin for the magnetization for the TRC with only even translations and only 
odd translations, respectively. In all cases except Zint,m,even, the dynamic exponent is a 
decreasing function of Lmin- For the magnetization, Zint,m appears to extrapolate to a value 
between 0.1 and 0.2 as Lmin oo while for the energy and staggered magnetization, the 
dynamic exponent appears to extrapolate to a value between 0.3 and 0.4. The small value of 
the dynamic exponent requires that odd translations and flips are included in the algorithm. 



From Fig. |TT] it is clear that the dynamic exponent is near 2 for the TRC algorithm with 
only even translations. 



Table |3 gives results of the weighted least squares fits for z for the smallest values of Lmin 
for which there is a reasonable confidence level. Since there is a general downward curvature 
in the log-log graphs, these numbers are likely to be overestimates of the asymptotic values. 
Thus, we can conclude that the asymptotic dynamic exponent for the TRC algorithm is 
likely to be less than 0.4 and is perhaps exactly zero. The dynamic exponent is apparently 
independent of the strength of the staggered field. For the case of the SW algorithm apphed 
to the two-dimensional Ising with no staggered field the best estimate is 2; = 0.25 ± 0.0100 
but the results are also consistent with logarithmic growth of relaxation times. The numbers 
for dynamic exponent for the SW appear to be smaller than for the TRC algorithm but this 
difference may simply reflect larger corrections to scaling in the case of the TRC . 
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V. DISCUSSION 



We studied the dynamics of the two-rephca cluster algorithm applied to the two- 
dimensional Ising model in a staggered field. We found that the dynamic exponent of 
the algorithm is either very small {z < 0.4) or zero ( r ~ logL) and that the dynamic 
exponent does not depend on the strength of the staggered field. A precise value of z could 
not be determined because of large corrections to scaling. We tested the importance of 
various optional components of the algorithm and found that an odd translation and global 
fiip of one replica relative to another is essential for achieving rapid equilibration. Without 
this component, z is near 2 so there is no qualitative improvement over the Metropolis al- 
gorithm. An odd translation and global fiip of one replica relative to the other allows for a 
large change of the total magnetization of the system with an acceptance fraction of 100%. 
Large changes in the global magnetization may also occur in the Swendsen-Wang algorithm 
in a field or via inactive fiips in the TRC algorithm but these fiips have a small acceptance 
fraction due to the staggered field. Unfortunately, the odd translation and fiip move is 
allowed because of a special symmetry of the Ising model in a staggered field. For more 
general Ising systems with translationally invariant fields, we expect performance similar to 
the TRC with even translations only. In this case, the autocorrelation time is significantly 
less than for the Metropolis algorithm but the dynamic exponent is about the same. While 
the two-replica approach is useful for these more general problems of Ising systems with 
fields, it does not constitute a method that overcomes critical slowing down except when 
additional symmetries are present that allow one replica to be fiipped relative to the other. 
Development of general methods for efficiently simulating critical spin systems with fields 
remains an open problem. 
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FIG. 1: Phase diagram of the two-dimensional staggered field Ising model, the three points on the 
critical line corresponds respectively to i7=0, 2, 4. 
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FIG. 2: Logarithm of the magnetization autocorrelation time Tint,m vs. logarithm of system size 
for -ff = 0, 2, 4. 
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FIG. 3: Magnetization autocorrelation time Tint,m vs. logarithm of system size L ioi H = 0, 2, 4. 
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FIG. 4: Logarithm of energy autocorrelation time Tint,£ vs. logarithm of system size L for H = 0, 
2, 4. 
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FIG. 5: Energy autocorrelation time Tint,e vs. logarithm of system size L ior H = 0, 2, 4. 




FIG. 6: Logarithm of staggered magnetization autocorrelation time Tint,s vs. logarithm of system 
size LioT H = 0, 2, 4. 
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FIG. 7: Staggered magnetization autocorrelation time Tint,s vs. logarithm of system size L for 
H = 0,2, 4. 
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FIG. 8: Finite size dynamic critical exponent for magnetization Zint,m vs. the reciprocal of the 



minimum size Lmm used in the fit. 
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FIG. 9: Finite size dynamic critical exponent for energy Zint,s vs. the reciprocal of the minimum 
size Lmin used in the fit. translations. 
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FIG. 10: Finite size dynamic critical exponent for staggered magnetization Zint,s vs. the reciprocal 
of the minimum size Lmm used in the fit. translations. 
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FIG. 11: Finite size dynamic critical exponent for magnetization Zint,m,even 

vs. the reciprocal of 
the minimum size Lmin used in the fit for the TRC with only even translations. 




FIG. 12: Finite size dynamic critical exponent for magnetization Zint,m,odd vs. the reciprocal of the 
minimum size Lmm used in the fit for the TRG with only odd translations. 
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TABLE I: Integrated autocorrelation times for the TRC algorithm for the magnetization of a single 
replica r^, the net staggered magnetization of both replicas and the energy of a single replica 
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H = 2 






H = A 




L(sizc) 




































' / III jn 




' i III jn 






' i III jn 




' iiii.S 


16 


10.7 ±0.1 


5.73 ± 0.09 


13.5 ±0.2 


10.8 ±0.2 


5.49 ± 0.08 


16.8 ± 0.2 


15.2 ±0.3 


5.7 ±0.1 


24 


17.0 ±0.2 


8.0 ±0.1 


20.2 ±0.3 


16.5 ±0.5 


8.0 ±0.2 


23.0 ±0.3 


22.4 ±0.4 


7.8 ±0.2 


32 


21.6 ±0.3 


9.8 ±0.1 


25.7 ±0.3 


21.8 ±0.3 


9.4 ±0.2 


28.1 ±0.4 


28.4 ± 0.6 


10.3 ±0.2 


40 


26.3 ± 0.4 


11.8 ±0.3 


29.2 ±0.5 


25.1 ±0.3 


11.4 ±0.2 


31.4 ±0.4 


31.7 ±0.8 


11.3 ±0.3 


48 


28.6 ±0.4 


12.9 ±0.3 


31.8 ±0.5 


28.4 ±0.8 


13.1 ±0.5 


34.4 ± 0.6 


35.6 ± 0.9 


12.3 ±0.3 


56 


30.6 ± 0.6 


14.7 ±0.5 


34.7 ±0.5 


28.7 ±0.6 


13.6 ±0.2 


37.5 ± 0.6 


39 ±1 


14.3 ± 0.4 


64 


34.2 ±0.6 


15.4 ±0.4 


36.2 ±0.6 


32 ± 1 


15.6 ±0.3 


38.1 ±0.5 


39.1 ±0.8 


15.3 ±0.6 


80 


37.3 ± 0.6 


16.2 ±0.4 


39.8 ±0.7 


37 ± 1 


16.7 ±0.5 


40.4 ± 0.8 


45 ±1 


16.7 ±0.8 


96 


39.5 ± 0.6 


19.0 ± 0.4 


42.6 ± 0.9 


40.2 ± 0.9 


18.0 ±0.4 


43 ± 1 


50 ±2 


17.9 ±0.3 


112 


41.7 ±0.8 


20.5 ± 0.7 


41.9 ±0.7 


43 ± 1 


19.0 ±0.7 


43.7 ±0.7 


50 ±1 


19.2 ±0.7 


128 


42.6 ±0.8 


20.4 ± 0.9 


43.7 ±0.8 


47 ± 1 


20.6 ±0.6 


45.8 ± 0.9 


54.4 ± 0.9 


19.4 ±0.6 


144 


44.6 ± 0.8 


21.5 ±0.7 


46.8 ±0.9 


47 ± 1 


20.8 ±0.4 


48 ± 1 


54 ± 1 


22.1 ±0.8 


160 


44.8 ± 0.7 


22.8 ±0.7 


44.7 ± 0.8 


55 ±2 


23.1 ±0.9 


48.5 ± 0.9 


60 ± 1 


22.3 ±0.6 


192 


46.2 ± 0.8 


24.0 ± 0.7 


47.2 ± 0.9 


56 ±2 


23 ±1 


50 ±1 


60 ±1 


25 ±1 


256 


50 ±1 


26.6 ± 0.9 


50.4 ±0.9 


58 ±2 


27.1 ±0.8 


52.2 ±0.9 


69 ±2 


1 26 ± 1.0 
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TABLE II: Integrated autocorrelation times for the magnetization of a single replica using 
TRC algorithm with odd translations only, Tint,m,odd and even translations only, Tint,m,even- 



L(sizej 


H- 


=0 


H- 


=2 


H 


=4 


'Tint,m,odd 


'^int,m,even 


'^int,m,odd 


'^int,m,even 


'^int,m,odd 


'^int,m,even 


16 


11.6 ±0.2 


12.6 ±0.1 


14.8 ± 0.4 


17.8 ±0.3 


19.0 ±0.3 


24.6 ±0.5 


24 


19.0 ±0.4 


25.7 ±0.3 


24.1 ±0.5 


34.9 ± 0.8 


28.4 ±0.5 


46 ± 1 


32 


25.9 ±0.3 


39.3 ±0.8 


30.3 ±0.4 


56 ±2 


35 ± 1 


77 ±5 


40 


32 ± 1 


58 ±2 


34.5 ± 0.8 


83 ±4 


44 ± 1 


109 ±4 


48 


36.6 ±1.0 


77 ±4 


39 ±1 


100 ±7 


44 ±2 


144 ±6 


56 


39.1 ±0.8 


106 ±9 


45 ± 1 


128 ±6 


49 ± 1 


188 ± 17 


64 


44 ± 1 


126 ±8 


48 ±1 


166 ± 11 


51 ±2 


203 ± 12 


80 


46 ± 1 


186 ± 18 


54 ±2 


283 ± 27 


55 ± 1 


434 ± 43 


96 


55 ±2 


— 


57 ± 1 


- 


59 ±2 


- 


112 


54 ±2 




60 ±3 




63 ±3 




128 


58 ±2 




62 ±3 




66 ±3 




144 


63 ±2 




65 ±2 




68 ±3 




160 


64 ± 1 




68 ±3 




70 ±3 




196 


68 ±3 




67 ±2 




75 ±4 




256 


69 ±2 




68 ±3 




77 ±4 
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TABLE III: Magnetization integrated autocorrelation times and CPU times for several algorithms 
for L = 80. 



A 1 P'orif.imn 


Integrated Autocorrelation Time 


CPU time 
(10~^ sec/sweep/spin) 


H = 


H = 2 


H = 4 


TRC 


37.3 ±0.6 


39.8 ± 0.7 


40 A ± 0.8 


3.1 


TRC 

odd translations only 


40 ± i 


04 ± Z 


r rr i 1 
00 ± i 


3.0 


TRC 

even translations only 


186 ± 18 


283 ± 27 


435 ± 43 


2.9 


TRC Sz inactive flips 
even translations only 


33.6 ±0.9 


246 ± 27 


372 ± 23 


4.6 


TRC 
no translations 


335 ± 18 


440 ± 24 


773 ± 47 


2.6 


Swendsen-Wang 


4.12 ±0.02 


4682 ± 173 


5707 ± 48 


1.3 


Metropolis 


928 ± 99 


1892 ± 158 


2959 ± 236 


1.1 
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TABLE IV: Estimated dynamic exponents together with minimum size used 
dence level for the TRC algorithm, TRC algorithm with odd translation only 
with even translation only. 



in the fit and confi- 
and TRC algorithm 



dynamic exponent z 


H = 


H = 2 


H = 4: 


Zint,m 


0.20 ± 0.03 

(112,86%) 


0.20 ± 0.02 

(80, 12%) 


0.23 ± 0.01 

(56, 75%) 


Zint,£ 


0.34 ± 0.04 
(80, 95%) 


0.40 ± 0.02 
(64,82%) 


0.42 ± 0.02 
(56, 60%) 


i,-^mini leVcV) 




0.42 ± 0.03 

(80,3%) 


0.34 ± 0.02 

(80, 14%) 


-^int,m,odd 


0.17 ±0.07 

(144,89%) 


0.23 ± 0.04 

(80, 60%) 


0.33 ± 0.02 

(40, 95%) 


Zint,m,even 


1.67 ±0.06 
(32, 96%) 


1.97 ±0.22 
(48, 71%) 
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